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ABSTRACT 


Equilibrium and kinetically limited performance, calculated using a one -dimensional 
model, are presented for methane and butene- 1 rocket fuels combined with fluorine- 
oxygen mixtures for chamber pressures of 100 and 300 psia (6. 89 5x10 5 and 2. 069X10 6 
N/m^) and an oxidant-fuel ratio range of 2. 5 to 6. Rate -controlling chemical reactions 
in the recombination process are identified and quantitatively rated in importance. Per- 
formance sensitivity to rate-constant uncertanity is studied in detail. It is shown that a 
one-dimensional kinetic flow model can be used to calculate performance in satisfactory 
agreement with the available experimental data and that no more than three reactions are 
rate controlling on the entire process. 


NONEQUILIBRIUM CALCULATIONS OF METHANE-FLUORINE - OXYGEN 
AND BUTENE- 1-FLUORINE - OXYGEN ROCKET PERFORMANCE 

by David A. Bittker 
Lewis Research Center 

SUMMARY 

Theoretical calculations of kinetically limited and equilibrium performance were 
carried out for the fuels methane and butene- 1 combined with fluorine-oxygen mixtures. 
The nonequilibrium results, obtained using a one-dimensional flow model, agree satis- 
factorily with the available experimental results for a chamber pressure of 100 psia 

c p 

(6. 895X10 3 N/m ) and a contoured nozzle having an area ratio of 40. The oxidant-fuel 
ratios used ranged from 2. 5 to 6 for oxidizers containing 54. 3 to 91.4 weight percent flu- 
orine. 

Two rating factors were used to identify the important chemical reactions occurring 
in the nozzle and to quantitatively rank them in importance as rate -controlling steps in the 
recombination process. Performance sensitivity to rate-constant uncertainty was studied 
in detail for the important reactions. Results showed that performance is controlled 
mainly by three reactions, even though several others must be included in the complete 
mechanism: H + H+ MsH 2 +M, H+F + M=HF + M, andH+OH+Ms HgO + M. 

The maximum kinetic performance for methane-fluorine - oxygen was 398 pound- 
seconds per pound (3903 (N)(sec)/kg) at 100 psia (6.895x10^ N/m^) chamber pressure 
using an oxidant containing 82. 6 weight percent fluorine. For butene- 1-fluorine - oxygen 
the maximum performance at this chamber pressure is 385 pound-seconds per pound 
(3776 (N)(sec)/kg) using an oxidant containing 70. 4 weight percent fluorine. 


INTRODUCTION 

The rocket performance characteristics of propellant systems containing fluorine 
have been under investigation in recent years by several laboratories. Theoretical and 
experimental performance of the hydrogen-fluorine system has been reported in refer- 
ences 1 to 3. Although the high performance level of the hydrogen-fluorine system 



makes it quite desirable for many missions, the practical problems of handling and 
storing liquid hydrogen have caused investigators to study the possibility of using less 
powerful, but more easily handled fuels combined with either fluorine, fluorine-oxygen 
mixtures, or the compound oxygen difluoride (OFg). The types of fuels considered are 
the low-molecular- weight hydrocarbons (e. g. , methane (CH^) , propane (CgH g ) or 
Butene- 1 (C 4 H g )), hydrazine (N 2 H 4 ), and the reactive compound diborane (B 2 H g ). Most of 
the possible combinations of the listed fuels and oxidizers are presently being considered 
for use as space-storable propellant systems. These are systems for which the propel- 
lants can be stored as liquids in outer space for long periods by using only insulating 
materials and heat shields. 

Figure 1 compares the performance of several space-storable systems with that of 
the best earth-storable and cryogenic liquid propellants in use today. This figure pre- 


System Chamber pressure, 

psia (N/m?) 

H 2 -F 2 60 (4. 137xl0 5 ) 

H 2 -0 2 300(2. 069xl0 6 ) 

B 2 H 6 -OF 2 100 (6. 895x10**) 

• B 2 H 6 -F 2 - 0 2 (70/30 mole %) 100 (6. 895xl0 5 ) 

N 2 H 4 -F 2 100(6.895xl0 5 ) 

CH 4 -F 2 - 0 2 (80/20 mole %) 100 (6.895xl0 5 ) 

CH 4 -OF 2 100(6.895xl0 5 ) 

Butene-1-F 2 - 0 2 100 (6. 895xl0 5 ) 

(66.7/33.3 mole %) 

Aerozine-50-N 2 0 4 100 (6. 895xl0 5 ) 



Equivalence ratio 


Figure 1. - Performance of rocket propellants. Equilibrium expansion to area ratio, 40. 
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sents equilibrium performance curves for six space-storable systems plus the hydrogen- 
fluorine, hydrogen-oxygen and Aerozine-50 - nitrogen tetroxide systems. Theoretical 
vacuum specific impulse is plotted as a function of equivalence ratio R defined by 



(Symbols are defined in appendix A. ) All calculations were performed using the method 
of Zeleznik and Gordon (ref. 4). It can be seen that the performance level of the space - 
storable propellants falls in the range of I = 380 to 435 pound-seconds per pound 

1 \ci C 

((lb) (sec) /lb) (3727 to 4266 (N)(sec)Ag) or midway between the limits of today’s highest 
and lowest energy systems. 

Two of the most promising systems under investigation at present are fluorine - 
oxygen mixtures combined with methane or butene- 1 as a fuel. Theoretical calculations 
for the methane combination were reported several years ago by Gordon and Kastner 
(ref. 5) . Recent work on these two systems was reported in reference 6, which gives a 
limited number of experimental data on altitude performance and which also presents 
theoretical performance curves based on the full equilibrium expansion assumption. Ref- 
erence 7 presents some theoretical performance calculations that consider the nonequi- 
librium (i. e. , finite reaction rate) nozzle expansion losses for the methane fluorine - 
oxygen system. The purpose of the current work is to present extensive nonequilibrium 
theoretical performance calculations for the methane -fluorine - oxygen and butene- 1- 
fluorine - oxygen systems. For each system, the results are compared with the limited 
experimental data available to determine how well the proposed mechanism and rates are 
able to predict the observed performance. The individual reactions in the proposed re- 
action mechanism are examined to determine which are most important to and controlling 
on the performance. The performance levels of both systerpfe are compared with each 
other so that their capabilities can be evaluated by engine designers. 

Conditions used for the methane -fluorine - oxygen system are chamber pressure, 

100 and 300 psia (6. 895X10^ and 2. 069x10® N/m 2 ); oxidant-fuel ratio O/F, 3 to 6; and 
weight percent fluorine in the oxidizer, 54. 3 to 91. 4. For the butene- 1-fluorine - oxygen 
system the chamber pressure is 100 psia (6. 895X10 5 N/m 2 ), O/F varies from 2. 5 to 5, 
and percent fluorine in the oxidizer ranges from 54. 3 to 82. 6 weight percent. 

^The specific impulse unit, (N)(sec)/kg, is identical to m/sec inasmuch as 1 N is 

9 

equal to 1 (kg)(m)/sec in the SI system. 
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PROCEDURE 


The complete investigation of rocket performance to be reported includes the fol- 
lowing steps: 


Equilibrium and Frozen Computations 

These computations were performed to obtain the limiting theoretical performance 
and to determine the significant species present for each system and, also, the equilib- 
rium starting conditions for the kinetic calculations. Computations showed that there 
was little or no effect on the results when certain trace species were purposely omitted 
from consideration. This fact allowed the selection of a set of significant species by 
using the criterion that omitting any species changes the equilibrium specific impulse by 
less than ±0. 5 (lb)(sec)/lb (±4. 9 (N)(sec)/kg) at an area ratio of 40. 

Kinetically Limited Computations 

With all significant species known, a reaction mechanism was then chosen for each of 
the systems under study. Many of the reactions were already known from previous in- 
vestigations; details of the selection procedure are given in the following section. A set 
of ’’reference" rate-constant values was chosen for all the postulated chemical reactions 
after a survey of the available literature. These rate constants were used to obtain ref- 
erence performance computations for all assigned conditions. The results of these cal- 
culations were analyzed to determine the most important reactions in each mechanism. 
Parametric rate variation studies were then made to find the effect of rate variations on 
performance. 

The nozzle used for all the kinetically limited computations is identical to the one 
used for the experimental measurements reported in reference 6. It is a truncated per- 
fect nozzle having a throat radius of 3 inches (7. 62 cm) and an area ratio of 40. Thrust 
level is approximately 5000 pounds (22 240 N) at 100 psia (6. 895x10^ N/m^) chamber 
pressure. 


Comparison with Experimental Results 

Experimental performance measurements for the methane - fluorine -oxygen system 
under altitude conditions are reported in reference 6. Table VI- 2 of that report tabulates 
uncorrected I yac and characteristic velocity efficiency r^*. In addition, the table lists 
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values of I corrected for heat transfer, friction, and nozzle divergence losses, as 
well as values of 77^* corrected for heat transfer and momentum losses. The I vac 
values were corrected for combustion inefficiency to obtain kinetically limited experi- 
mental specific impulse data that were compared with the computed theoretical results. 
The equation used is 


j _ ^vac^corr ^2) 

^C^corr 

THEORY AND INPUT DATA 
Method for Kinetic Calculations 

The problem of calculating kinetically limited rocket performance is part of the more 
general problem of calculating the progress of a complex set of chemical reactions in a 
flow system. One must solve the usual differential equations of fluid motion (conser- 
vation of mass, momentum, and energy) simultaneously with the equations for the 
changes of concentration due to chemical reaction. The particular system of equations 
used in the present work assumes one-dimensional flow of uniformly mixed, thermally 
perfect, reacting gases. The reactants are assumed to burn with 100 percent combustion 
efficiency in the chamber and all flow losses due to the effects of viscosity and other 
transport phenomena are neglected. The specific equations and computer program used 
for the computations were developed under NASA contract and are described in refer- 
ence 8. They are the same as those used for the hydrogen-fluorine calculations de- 
scribed in reference 2. To perform these kinetic computations a set of reactions that 
accurately describe the chemical changes occurring in the flow system must be 
written. All reactions are considered reversible, and the law of microscopic reversi- 
bility is used to obtain the reverse rate constant for any reaction from its forward rate 
constant and equilibrium constant. 


Reaction Mechanism and Species for Hydrocarbon-Fluorine - Oxygen System 

Early theoretical work by Gordon and Wilkins (ref. 9) showed that the carbon atom 
combines preferentially with the oxygen atom and hydrogen with fluorine in a 
hydrocarbon-fluorine - oxygen combustion. The two major products of combustion are 
carbon monoxide (CO) and hydrogen fluoride (HF) . Reference 9 also showed that best 
theoretical performance was obtained when the number of carbon atoms is the same as 


5 


the number of oxygen atoms in the mixture and the number of hydrogen atoms equals the 
number of fluorine atoms. This means that the carbon to hydrogen atom ratio of the fuel 
will be the same as the oxygen to fluorine atom ratio in the oxidant. For example, the 
following reaction can be written for the greatest performance of the methane -fluorine - 
oxygen system: 


CH 4 + i 0 2 + 2F 2 - CO + 4HF 


(A) 


This reaction has an oxidant containing 82. 6 weight percent fluorine and an oxidant-fuel 
ratio of 5. 75. For any other assigned oxidant composition the best performance will be 
obtained for one particular O/F value. A simple approximate method of calculating this 
value is also given in reference 9, and the equations are presented for reference in ap- 
pendix B. Results of the present equilibrium calculations showed that two different sets 
of significant species are obtained as a function of the percent fluorine in the oxidant and 
the assigned O/F value. These variables control the atom ratios in the mixture. By 
using the oxygen to carbon atom ratio a defined in appendix B, the two species lists and 
corresponding chemical reaction mechanisms can be classified as follows: 

Mechanism A: oxygen-to-carbon-atom ratio less than 1 (low O/F) and high fluorine- 
to-oxygen ratio. - For this situation CO is the only oxygen containing species, and the 
excess carbon atoms can combine with fluorine and hydrogen atoms. The complete list 
of significant species needed to satisfy the selection criterion described eariler is gase- 
ous carbon (C), hydrogen atom (H), fluorine atom (F), carbon monoxide (CO), hydrogen 
fluoride (HF), difluoroacetylene (CgFg), acetylene (CgHg) and methylidyne radical (CH). 
No solid carbon or other condensed products were found. The reaction mechanism as- 
sumed for these species considers both CO and CgF 2 as inerts and is as follows: 


H + F+ M £ HF + M 
AHggs = -136 kcal/mole f 
= -569 kj/mole 


H + H+Mf H 2 +M 
AHggs = - 104 kcal/mole 
= -435 kj/mole 


> 


j 


(I) 


(II) 
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* + H 2 *° 


HF+ H 


AHggs = -32 kcal/mole 


= -134 kj/mole 


k 4 

CH + CH 5? C 2 H 2 


A 


AH ggg = -231 kcal/mole 
= -964 kj/mole 
kc 

C + H + M ? 3 CH + M 


A 


AHggg = -81 kcal/mole 


= -339 kj/mole 


(m) 


(IV) 


(V) 


Me ch anism B: oxygen- carbon ratio greater than one (hig h O/F) and low fluorine to 
oxygen ratio . - In this case there is excess oxygen present to combine with hydrogen atom 
as well as carbon. The fluorine forms only HF and no condensed species were found. The 
complete list of species is H, O, F, CO, C0 2 , H 2 , HF, 0 2 , OH, and HgO. The reaction 
mechanism for these species consists of reactions I, II, and in plus the following: 


CO + O + M £ C0 2 + M 
AHggg = -127 kcal/mole 
= -531 kJ/mole 


> 


CO + OH s' C0 2 + H 
AH 29 g = -25 kcal/mole 
= -105 kj/mole 


(VI) 


(vn) 


7 



AHggs = -119 kcal/mole ). 


= -498 kj/mole 


(vni) 


o + o+ m^o 2 +m 


AH 298 = ~ 119 kcal / mole y 


= -498 kj/mole 
O+H+M *° OH + M 


AH 298 = “ 192 kcal / m °l e r 


= -427 kj/mole 


(IX) 


(X) 


H 2 + OH b 1 H 2 0 + H 


AHggg = -15 kcal/mole 


= -63 kj/mole 


(XI) 


H+ Or 


L 12 


OH + O 


AH 298 = +17 kca l/ m °l e f 


= +71 kJ/mole 




K 1S 

O + H 2 ; OH + H 

AH 298 = +2 kca l/ m °l e f 


(XII) 


(xm) 


= +8 kJ/mole 



Experimental data are available for most of the chemical reactions written. After a 
survey of the current literature presenting both theoretical and experimental rate- 
constant data, the rate-constant expressions shown in table I were chosen as the set of 
values to be used. These rate constants will be referred to as the reference rates. It 
should be pointed out that, although the rate constants for reactions I, n, and HI are 
taken from reference 2, they are not necessarily the ''assigned’' rates used in that work. 
In fact, only k 9 in table I is the same as the assigned rate constant of reference 2. The 


TABLE I. - RATE CONSTANTS FOR ME THANE -FLUORINE 
OXYGEN SYSTEM 3 -’ b 


Reaction 


Reaction Pre exponential 
number factor, 

A 


H+F+IVfctHF+M 

I 

H+H+ M^H 2 + M 

II 

f+h 2 ^hf+h 

m 

CH+CH=tC 2 H 2 

rv 

C+H+M^CH+M 

V 

ccm-o+m^co 2 +m 

VI 

CCM-OH^C0 2 +H 

VII 

h+oh+m^h 2 o+m 

VIII 

O+O+M^Og-i-M 

IX 

OfH+M^OH+M 

X 

H 2 +OH^H 2 Ot-H 

XI 

h+o 2 ^oh+o 

XII 

Oi-H 2 ^OH+H 

xm 

f+oh^hf+o 

XIV 


3.0X10 
7. 5X10 

2. 04X10 
5X10 
2X10 
7.0X10 

3. 3X10 

4. 5X10 
8. 15X10 

4.0X10 
2. 3X10 
2. 2x10 
1. 2X10 
2.9X10 


19 

18 

13 

11 

16 

16 

11 

21 

18 

18 

13 

14 
13 
12 


Temperature 

Activation 

Reference 

exponent, 

n 

energy, 

E a’ 

cal/mole 


-1.0 

0 

2 

-1.0 

0 

2 

-. 12 

3 750.0 

2 

.5 

6 000. 0 

16 

-.5 

0 

16 

-1.0 

0 

C 15 

0 

600.0 

17 

-1.5 

0 

18 

-1.22 

0 

19 

-1.0 

0 

20 

0 

5 200.0 

17 

0 

16 600. 0 

21 

0 

9 200.0 

22 

.68 

200. 0 

11 


a Rate constant expression is k = AT n exp(-E a /RT). 

b Units of k are cm 3 /(mole) (sec) for two-body reactions and cm 6 /(mole 2 )(sec) 


for three-body reactions. 

c Data from ref. 15 fitted to the equation k = AT"*. 


expressions for kj and k g given in table I are derived from recent experimental re- 
sults (ref. 10) and are used in the present work as reference rate constants. 


RESULTS AND DISCUSSION 
Methane-Fluorine - Oxygen System 

Performance results. - Kinetic performance has been calculated for a wide range of 

— — — — — 5 g 

operating conditions. At chamber pressures of 100 and 300 psia(6. 895x10 and 2. 069x10 
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Vacuum specific impulse, I vac , (N)(sec)/kg 


3800 - 
3600 - 
3400 - 
3200 L- 

4000 r~ 
3800 — 
3600 - 
3400 L 

j 

4200 r- 

! 

4000 - j 
3800 — 
3600 — 
3400 — 

4000 1 — 
3800 — 
3600 — 
3400 — 



Chamber pressure, 
psia (N/m 2 ) 

300.(2. 069xl0 6 ) 
100.(6. 895X10 5 ) 

300(2.069xl0 6 ) _ i 
I I I 
100 (6.895x10 s ) 


300 (2. 069xl0 6 ) — 
100 (6.895x10 s )- 

I I I 


(a) Fluorine content in oxidant, 54. 3 weight percent. Stoichiometric O/F, 5. 82. 




(c) Fluorine content in oxidant, 82. 6 weight percent; stoichiometric O/F, 7.65. 


300 (2.069xl0 6 )~ 
100 (6.895xl0 5 )- 
300 ( 2. 069xl0 6 )_ 

l l i 

100 (6.895xl0 5 )- 


300 (2. 069xl0 6 ) 
100 16. 895x10 s ) 

6.2 


(d) Fluorine content in oxidant, 91.4 weight percent; stoichiometric O/F, 8.48. 
Figure 2. - Methane-fluorine - oxygen performance. Area ratio, 40. 
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N/ m 2 ) the oxidant-fuel ratio (O/F) was varied from 3 to 6 for four oxidants, containing 
50, 70, 80, and 90 mole percent fluorine. Weight percents of fluorine in these oxidant 
mixtures were 54. 3, 73.5, 82.6, and 91.4, respectively. 

A complete summary of equilibrium, kinetic, and frozen performance is shown in 
figure 2. In these graphs the indicated stoichiometric O/F values are those for conver- 
sion to the final products HF, CC> 2 , and HgO. The highest level of performance is shown 
with the oxidant containing 82. 6 weight percent fluorine. Figure 2(c) shows that the max- 
imum kinetic specific impulse is 398 (lb) (sec) /lb (3903 (N)(sec) /kg) at O/F = 5. 1 for 
100 psia (6. 895X10 5 N/m 2 ) chamber pressure and 412 (lb)(sec)/lb(4040 (N)(sec)/kg) at 
O/F = 5. 6 for 300 psia (2. 069x10^ N/m 2 ). To better visualize the effect of oxidant com- 
position, performance has been plotted against percent fluorine in the oxidant for several 
O/F values in figure 3. Both the equilibrium results (fig. 3(a)) and the kinetic results 


Oxidant-fuel 

ratio, 

O/F 

6 

5. 75 



Figure 3. - Effect of oxidant composition on methane-fluorine - oxygen performance. 
Chamber pressure, 100 psia (6,895xl0 -5 N/m 2 U final area ratio, 40. 



(fig. 3(b)) show that highest performance is obtained for oxidants containing 82 to 86 per- 
cent fluorine. This fact agrees with reaction (A), written previously, and confirms the 
conclusion made from figure 2. It should be noticed that figure 2(c) shows maximum ki- 
netic performance occurring at an O/F very close to 5 for low chamber pressure. This 
value is much lower than the O/F of 5. 75 obtained from reaction (A). At the high 
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TABLE II. - IMPULSE RECOVERY PARAMETERS 


FOR ME THANE -OXYGEN-FLUORINE 
[Oxidant-fuel ratio, 5. ] 


(a) U. S. Customary units 


Fluorine content 
in oxidant, 
wt. % 

Chamber pressure, 

P c , psia 

100 

300 

100 

300 

Impulse recovery factor, 
F 

(a) 

Impulse loss, 
I eq “ 

(Ib)(sec)/Ib 

54.3 

0.660 

0.817 

18.2 

9. 1 

73.5 

.665 

.852 

18.0 

7. 1 

82.6 

.694 

.850 

15.5 

6.7 

91.4 

.695 

.838 

14.8 

7.0 


(b) SI units 


Fluorine content 
in oxidant, 
wt. % 

o 

Chamber pressure, P c , N/m 

6. 895X10® 

2. 069X10 6 

6.895X10 5 

2. 069x10® 

Impulse recovery factor, 
F 

Impulse loss, 

! eq ‘ I R’ 
(N) (sec) /kg 

54.3 

0.660 

0.817 

178 

89 

73.5 

.665 

.852 

177 

70 

82.6 

.694 

.850 

152 

66 

91.4 

.695 

.838 

145 

69 


^^Ik-IfP/deq-y- 


chamber pressure, however, maximum performance occurs close to O/F = 5. 75. One 
explanation of this behavior is the deviation from the assumption that the nozzle gas con- 
sists of only the species HF and CO. This assumption is good for an equilibrium condi- 
tion but is much less valid for a nonequilibrium situation with much slower recombination 
rates for the formation of HF. The increase of these recombination rates with increased 
pressure explains why the kinetic performance curve for 300 psia (2. 069X10® N/m 2 ) 
chamber pressure has its maximum at an O/F value closer to 5. 75 than does the low- 
pressure curve. From the practical viewpoint, achieving maximum performance at a 
lower O/F means that less oxidant has to be carried along to perform a given mission. 

The curves of figure 2 show that increasing the chamber pressure from 100 to 300 
psia (6. 895X10 5 to 2. 069X10® N/m 2 ) has a small effect on equilibrium performance but 
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very significantly increases the reaction rate limited performance in all cases. This is 
because the largest effect of the pressure increase is to increase the rate of three-body 
collisions, and, thus, the overall rates of three-body recombination reactions. These 
are the reactions that release most of the energy into the expanding gases and are there- 
fore rate controlling for kinetically limited performance. As a measure of performance 
losses due to finite reaction rates, we will use the impulse recovery parameter, defined 

by 


F = 



( 3 ) 


For comparison purposes the values of F are calculated for all the conditions in fig- 
ure 2 at an O/F of 5. Results are given in table n and show quantitatively the percent- 



Figure 4. - Theoretical performance against nozzle area ratio for 
methane-fluorine - oxygen. Chamber pressure, 100 psia 
(6. 895x10* oxidant-fuel ratio, 5.75; fluorine content in ox- 
idizer, 82. 6 weight percent. 
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age decrease of kinetic losses with increase in chamber pressure as well as the absolute 

performance loss I Q _ - L. On an absolute basis, at O/F = 5, the kinetic loss I - I, 
eq k 5 2 eq k 

at a chamber pressure of 100 psia (6. 895x10 N/m ) is about twice the loss at 300 psia 

(2.069X10 6 N/m 2 ). 

It is also of interest to see how the kinetic performance loss depends on position in 
the nozzle. A typical plot of the three types of theoretical performance as a function of 
supersonic nozzle area ratio is shown in figure 4. The 82. 6 percent oxidizer is used at 
an O/F of 5. 75 and a chamber pressure of 100 psia (6. 895X10 5 N/m 2 ). At an area ratio 
of 2, the kinetic loss I - I k is 1 impulse unit or less than 1/2 percent of the equilib- 
rium value. A 1-percent loss in performance occurs at an area ratio of 3. 5. For an 
area ratio of 40 the kinetic loss is about 20 impulse units or 5 percent of the equilibrium 
value. 

Effect of reaction rate variation. - It is possible to obtain a qualitative idea about 
which reactions may be the most important in determining the calculated specific im- 
pulse by examining the changes in various species concentrations downstream of the 
rocket nozzle throat. Some typical examples of these concentration changes are shown in 
table HI. Species mole fractions near the throat of the nozzle and at an expansion ratio 
of 40 are tabulated at O/F = 3 and 5. 75 for two oxidizers. 

For both oxidizers at O/F = 3, the only significant changes are a decrease in hydro- 
gen atom concentration and an increase in molecular hydrogen concentration. At 
O/F = 5. 75 both hydrogen and fluorine atom concentrations decrease, whereas HF in- 
creases. The molecular hydrogen concentration now decreases instead of increasing. 
From this information one might expect that the hydrogen atom recombination is most 
rate controlling at O/F = 3, but several possible reactions, for example, (I), (n), (ID), 
and (XI), could be rate controlling at O/F = 5. 75. A more systematic approach to 
measuring the importance of each chemical reaction has been suggested in reference 11, 
which assumes the importance of a reaction to be proportional to the net volumetric con- 
version rate, that is, the extent to which the reaction proceeds. The general formula for 
this quantity can be found in texts such as the work of Penner (ref. 12). The rating factor 
to be used is given by the formula 


X 

v v 


N 

k f IT 

,3 = ,] i=l 


c i 1] “ k b i 

1 b ’J i=l 


vU 


ij 


(4) 
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TABLE in. - COMPOSITION OF EXHAUST GAS FOR 


ME THANE -FLUORINE - OXYGEN PROPELLANT 
[Chamber pressure, 100 psia (6. 895x10® N/m^). ] 
(a) Fluorine content in oxidizer, 82. 6 weight precent 


Species 


Oxidant fuel ratio 




3 

5.75 


Throat 

Area ratio, 40 

Throat 

Area ratio, 40 



Species mole fraction 

C 

0.00009 

0.00009 

0.00011 

0.00011 

CO 

. 14158 

. 14339 

. 17971 

. 18687 

C 2 F 2 

.05354 

.05422 

.00105 

.00109 

CH 

. 00003 

. 00004 

.000001 

. 000003 

c 2 h 2 

. 01093 

.01107 

5.6X10' 9 

.000001 

.000001 

F 

. 00096 

. 10380 

.05527 

H 

. 06673 

.04307 

. 08858 

.06625 

h 2 

. 26783 

. 28301 

.01294 

. 00005 

HF 

.45830 

.46512 

.61364 

.69035 

(b) Fluorine content in oxidizer, 73. 5 weight percent 


0.00001 

0.00001 






CO 

co 2 

. 20796 

. 20998 

0. 17503 
. 00927 

0. 16749 
.02335 

.02231 

. 02253 




C 2 F 2 

CH 



.00001 

.00001 



C 2 H 2 

F 

. 00405 
. 00070 

. 00409 
1.0X10" 10 

.05210 

. 00449 

H 

.05530 

.03715 

.07381 

. 07773 

h 2 

. 26972 

. 28133 

.01724 

.00345 

HF 

H 2° 

.43993 

.44491 

.59595 
. 00760 

. 66654 
. 00584 



O 

. 04365 

.03416 




. 00756 

. 01644 



OH 

.01779 

.00051 





TABLE IV. - REACTION RATING PARAMETERS FOR 


METHANE-FLUORINE - OXYGEN FUEL 
[Chamber pressure, 100 psia (6. 895x10^ N/m^). ] 


(a) Fluorine content in oxidizer, 82. 6 weight percent 


Reaction 

number 

Reaction 

Mixture ratio 

3 

6 

Reaction rating parameter 

x j 

Rank 

X H,j 

Rank 

x i 

Rank 

X H, j 

Rank 

(ID 

H+H+M = H 2 +M 

958.0 

1 

-99 610 

1 

239 

3 

-24 850 

2 

(i) 

H+F+M = HF+M 

56. 1 

2 

-7627 

2 

1710 

1 

-232 600 

1 

(in) 

F+H 2 = HF+H 

22.9 

3 

-731 

3 

573 

2 

-18 350 

3 

(V) 

C+H+M a CH+M 

0. 168 

4 

-1.36 

4 


-- 


-- 

(IV) 

CH+CH = C 2 H 2 

8.9X10" 6 

5 

-2. 06X10 -3 

5 


— 


-- 

(vra) 

H+OH+M = H 2 Of M 


— 


— 

33. 2 

5 

-3946 

4 

(X) 

O+H+M = OH+M 


-- 


-- 

9.85 

7 

-1004 

5 

(XI) 

h 2 +oh = h 2 o+h 


-- 


— 

-34. 1 

4 

+511 

6 

(VII) 

COfOH = C0 2 +H 


— 


-- 

14. 1 

6 

-353 

7 

(VD 

COfO+M = C0 2 +M 


— 


— 

0. 579 

10 

-74 

8 

(EX) 

0+0+ M = 0 2 +M 


— 


— 

0.258 

11 

-31 

9 

(XII) 

H+0 2 = OH+O 


— 


— 

-1.38 

9 

-24 

10 

(xni) 

0+H 2 = OH+H 


-- 


— 

-9. 14 

8 

00 

t-H 

l 

11 


(b) Fluorine content in oxidizer, 73. 5 weight percent 


(id 

H+H+M = H 2 +M 

198 

6 

-20 590 

4 





a) 

H+F+M = HF+M 

916 

1 

-124 600 

1 





(m) 

F+H 2 = HF+H 

769 

2 

-24 600 

3 





(vm) 

H+OH+M = H 2 0+M 

476 

3 

-56 640 

2 





(X) 

O+H+M = OH+M 

89. 5 

9 

-9 129 

5 





(XI) 

h 2 +oh = h 2 o+h 

-445 

4 

+6 675 

6 





(VII) 

CO+OH = C0 2 +H 

200 

5 

-5 000 

7 





(VI) 

CO+O+M = C0 2 +M 

5.5 

11 

-699 

10 





(K) 

0+0+ M = 0 2 +M 

24.8 

10 

-2 951 

8 





(XII) 

H+0 2 = OH+O 

-142 

7 

-2 414 

9 





(xm) 

0+H 2 = OH+H 

97 

8 

+ 194 

11 








In this equation X - is the volumetric conversion rate tot the reaction 

V, j 


N k N 

V f » j V , 

V 1 k b,j La V 1 

i=l i=l 


in 

where IVL designates the i in species in the reacting system. The sign of indicates 
the direction in which the reaction is actually proceeding. A positive X^ means the re- 
action is actually producing products from reactants as written. 

Performance, however, depends not only on a net reaction rate, but also on the 
amount of heat released into or absorbed from the system by each reaction. In the pres- 
ent work, therefore, a more meaningful parameter has also been used, namely, 


X H,j = X j< iH 298> 


( 5 ) 


where AH 


298 


is the molar heat of reaction at 298 K for the forward reaction as written. 


The units of AHggg are arbitrarily taken as kilocalories per mole. If the usual sign con- 
vention is used for AHggg (A negative sign means heat is released by the reaction. ), then 
the sign of X H,j will always be negative if heat is released into the gas and positive if 
heat is absorbed from the gas. For example, if a reaction, as written, is exothermic, 
but actually proceeds in the reverse direction, Xj is negative and AHggg is negative. 
Then X H, j is positive, indicating that this reaction absorbs heat from the gas; there- 
fore, X H • gives both the magnitude and direction of the enthalpy change rate caused by 
a particular reaction. 

Values for both rating parameters have been calculated for the reactions in the 
methane -fluorine - oxygen expansion for three conditions, and are given in table IV. 

These show that either X. or X™ . selects reactions (I), (II), and either (HI) or (VUI) 
as the most important at all operating conditions. To determine the actual sensitivity of 
performance to uncertainty in the rate-constant values for reactions (I), (H), and (HI), 
rate variation calculations were performed in which each of the three reaction rates was 
individually changed. The rate variation factors were chosen to account for experimen- 
tal uncertainties and third body efficiency effects. The effect of reaction II (the H atom re- 
combination) on performance is shown in figure 5 for two oxidant compositions, 73. 5 and 
82. 6 weight percent fluorine. The rate constant variation is from five times the refer- 
ence value to one fourth of that value. It can be seen that either increasing or decreasing 
kg has a significant effect on performance over almost the entire range of O/F values 
for both oxidants. Even at O/F = 6, where the effect is smallest, increasing kg causes 
almost a 1 percent increase in specific impulse for both oxidants. 
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Relative rate of H atom 
recombination 
reaction 

® (equilibrium) 

— 0.25 

0 (frozen) 


O Corrected experi- 
mental data 



(b) Fluorine content in oxidant, 82.6 weight percent. 

Figure 5. - Performance of methane-fluorine - oxygen: effect of 
hydrogen atom recombination rate. Chamber pressure, 100 psia 
(6.895X10 5 N/m2) ; area ratio, 40. 


The effect of changing only kj, the rate constant for the H + F recombination re- 
action, is shown in figure 6 for the same conditions used in figure 5. For this reaction, 
kj was increased to four times the reference value and decreased to one-fourth of this 
value. The results in figure 6 show that changing k^ has negligible effect on perform- 
ance at O/F values less than 4, but it does have a significant effect at the higher O/F 
values. The maximum change is about 8 impulse units or 2 percent for the 82. 6 weight 
percent fluorine oxidizer. 
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Relative rate of HF 
atom recombination 
reaction 




> 


°> (equilibrium) 

0.25 

0 (frozen) 


O Corrected experi- 
mental data 



(a) Fluorine content in oxidant, 73. 5 weight percent. 





Oxidant-fuel ratio, O/F 


(b) Fluorine content in oxidant, 82.6 weight percent. 


Figure 6. - Performance of methane-fluorine - oxygen: effect of 
hydrogen atom plus fluorine atom recombination rate. Chamber 
pressure, 100 psia (6.895X10 5 N/m 2 ); area ratio, 40. 


Table V shows the effect of changing kg for the reaction F + Hg? HF + H. The 
value of kg was changed by factors of 10 and 0.1. The data in table V show that per- 
formance is much less sensitive to variation in the rate of reaction (III) than to variation 
of the rates of reactions (I) and (II) . Neither of the rate changes causes more than a 
0. 6 unit change in specific impulse. Any change which is no greater than 0. 5 impulse 
unit in magnitude was considered negligible in this work. From the last column in 
table V it is evident that removing reaction (HI) entirely from the mechanism at low O/F 
values changes performance negligibly. Its removal has a quite significant effect at high 
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TABLE V. - EFFECT ON PERFORMANCE OF RATE 




OF THE REACTION F+H 2 s' 3 HF+H 


Fluorine content 

Oxidant -fuel 

Rate constant, 

kg, equal to - 

in oxidant, 
wt. % 

ratio, 

O/F 

k r 

5k R 

0.1 k R 

0 



Vacuum specific impulse, I , (lb) (sec) /lb 

VdC 

73.5 

3.5 

377.7 

377. 7 

377.7 

377.6 

73.5 

6 

384.3 

384.9 

384.2 

383.5 

82.6 

4.0 

385.0 

385.0 

385.0 

384.8 

82.6 

5.667 

397. 2 

397.2 

396.8 

394.7 


a kp = Reference rate constant. 


O/F, however, where the concentration of the fluorine atom is quite high. The matter of 
removing a reaction will be discussed further at the end of this section. It should be 
noticed from table IV that j correctly predicts that performance should be more de- 
pendent on the rate of the H + H+ MjHj+M reaction at high O/F than on the rate of 
F + Hg s HF + H, whereas does not. 

Table IV(b) shows that reaction (Vin), H + OH + M z HgO + M, should be important 
at high O/F for the 73. 5 -percent oxidizer. Inasmuch as this reaction is stated to have 
a rather large uncertainty in reference 11, the effect of its rate on performance was in- 
vestigated for three conditions of high oxidant-fuel ratio. The results, as shown in 
table VI, agree with the information in table IV. Performance is sensitive to the rate 
constant of reaction (Vin) for the 73. 5 weight percent oxidizer, but insensitive to kg for 
the 82. 6 weight percent oxidizer. The value kg = 0. 10 k^ is the rate constant suggested 
in reference 11 and is identical to the rate constant for H + OH + Arj HgO + Ar inferred 


TABLE VI. - EFFECT ON PERFORMANCE OF RATE OF 


THE REACTION H+OH+M H 2 CM-M 


Fluorine content 

Oxidant-fuel 

Rate constant, 

kg, equal to - 

in oxidant, 

ratio, 

k R 




wt. % 

O/F 

10 kR 

0.1 k R 

0 



Vacuum specific impulse, I yac , (lb) (sec) /kg 

73.5 

5.667 

' 

387.8 

397. 5 

385.3 

384.9 

73.5 

6 

384.3 

393.4 

382. 5 

382.2 

82.6 

6 

394. 1 

394.9 

393.9 

393.8 


a kp = Reference rate constant. 
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from the experimental shock tube data reported in references 13 and 14. Our uncertainty 
limits for this reaction are considerably narrower than those estimated in reference 11. 
The data in table VI show that reaction (VUE) would, because of its rate-constant uncer- 
tainty, be an important reaction to consider, if an engine were operated at high O/F 
values with the 73. 5 weight percent fluorine oxidant. 

Rate variation calculations are also presented for reaction (VI) (CO + 0+ Mr 
C0 2 + M), which has very small Xj and Xg . factors but a large degree of uncertainty 
in its rate constant. The reference rate expression in table I was obtained by curve fit— 

* ting the data of reference 15. Its numerical values are much lower, for a given tempera- 
ture, than those given by the expression in reference 11. Consequently, rate variation 

• factors of 200 and 0. 5 were used to obtain the results shown in table VII. The data indi- 

TABLE vn. - EFFECT ON PERFORMANCE OF VARIATION 

k fi 

OF RATE OF THE REACTION COfCM-M J 3 CO a +M 


[Oxidant fuel ratio, 6. ] 


Fluorine content 

Rate constant, 

kg, equal to - 

in oxidant, 
wt. % 

k R 

200 k R 

0.5k R 

0 


Vacuum specific impulse, , (lb) (sec) /lb 

73.5 

384.3 

387.9 

384.7 

384.7 

82.6 

394. 1 

394.0 

394.0 

394.0 


a k R = reference rate constant. 


cate that performance is not sensitive to the value of kg when the 82. 6 weight percent 
fluorine oxidant is used. Performance is affected by kg, however, for the 73. 5 weight 
percent fluorine oxidant. In comparing the kg expression in table I with the one given 
in reference 11, it was concluded that values from the table I may be slightly too low, 
because of third body efficiency effects, but that the reference 11 expression probably 
gives values that are much too high. In the high-temperature shock tube work of refer- 
ence 15, care was taken to remove all traces of water, which is known to act as a cata* 
lyst and could be the cause for erroneously high rate constants for the CO + O recombi- 
nation process previously reported. 

Effect of reaction F + OH ^ HF + O . - The work of reference 11 indicates that one 
reaction, which has not been considered herein, may be important and should be included 
in the reaction scheme for methane -fluorine - oxygen. This is the reaction 
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F + OH ? HF + O 
Hggs = -33 kcal/mole > 
= -138 kJ/mole 


(XIV) 


Calculations were made for two operating conditions at which this reaction should show 
its maximum effect. These are for the 73. 5 and 82. 6 weight percent fluorine oxidizer at 

C n 

O/F = 6 and a chamber pressure of 100 psia (6. 895x10° N/ni ) . Results are shown in 
table VHI. Because reaction (VI) affects the oxygen atom concentration, the computations 
were done for two values of kg, namely, the reference value of table I and 200 times this 


TABLE vm. - PERFORMANCE CHANGE DUE TO 
THE REACTION F+OH k | 4 HF+O 


[Oxidant-fuel ratio, 6; chamber pressure, 100 psia 
(6.895X10 5 N/m 2 ).] 


Rate 

Fluorine content 

Vacuum specific impulse change, 

Case 

constant 

in oxidizer, 

at b 

^vac’ 


for 

wt. % 

(lb)(sec)/lb 


reaction 




(VI) 




\ 

k R 

73.5 

+0.23 

1 

82.6 

+ .03 

2 

200 k R 

73.5 

+0.51 

3 

82.6 

+ .05 

4 


a Reference rate of constant. 

k^vac = Vac reaction (XIV)) - I vac (without reaction (XIV)). 


value. The conclusion from these results is that the addition of reaction XIV has negli- 
gible effect on performance for the two conditions investigated. The largest increase in 
I caused by this reaction is 0. 51 impulse unit, when the high rate for the 

V3.C 

CO + 0+ Mj COg + M reaction is used. The X H, j rating parameter values for cases 1 
and 3 of table VHI are shown in table IX. It is seen that reaction (XIV) ranks last in im- 
portance if the reference value is used for kg and sixth out of 12 if the high kg value 
is used. Thus, even when reaction (XIV) ranks sixth out of 12, its addition to the mech- 
anism changes specific impulse only by the assumed limit for an unimportant reaction. 
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TABLE IX. - REACTION RANKING PARAMETERS WITH THE 


REACTION F+OHsHF+O INCLUDED IN MECHANISM 


[Chamber pressure, 100 psia (6. 895x10’’ N/in^); oxidant-fuel ratio, 
6; fluorine content in oxidant, 73. 5 wt. %. ] 


Reaction 

Reaction 

Case 1: = 

k R 

Case 3: k fi = 200 k R 

number 





Reaction rating 

Rank 

Reaction rating 

Rank 



parameter, 


parameter. 




X H,j 


X H,j 


I 

H+F+MzHF+M 

-121 200 

1 

-87 970 

1 

vm 

h+oh+m^o+m 

-56 450 

2 

-40 920 

3 

n 

h+h+m=h 2 +m 

-20 570 

3 

-14 840 

4 

m 

f+h 2 ^hf+h 

-17 330 

4 

-13 400 

5 

X 

O+H+M^OH+M 

-9 584 

5 

-6 942 

8 

XI 

h 2 +oh=h 2 cm-h 

+6 771 

6 

+4 744 

9 

vn 

cooh^co 2 +h 

-4 723 

7 

+7 667 

7 

IX 

0+0+Mi0 2 +M 

-3 101 

8 

-2 248 

11 

xn 

h+o 2 ^oh+o 

-2 332 

9 

-2 292 

10 

VI 

C0+0+M^tC0 2 +M 

-726 

10 

-73 850 

2 

xin 

o+h 2 =oh+h 

+660 

11 

+508 

12 

XIV 

F+0H=KF+0 

L.. 

-238 

12 

-12 070 

6 


In summary, rate variation calculations have shown that performance of the methane- 
fluorine - oxygen system is mainly controlled by the hydrogen atom recombination, the 
H + F recombination, and, under some conditions, the H + OH recombination. Next in 
importance is the reaction F + Hg ^ HF + H. Although this reaction must usually be 
present, variation of its rate constant by one order of magnitude causes negligible change 
in computed specific impulse. This last result is consistent with the fact that bimolecu- 
lar reactions are seldom rate controlling in the recombination process because of their 
high rates and relatively low heats of reaction. They generate the atomic and free radi- 
cal species that participate in the rate -controlling three-body recombinations. 

One more matter remains to be discussed. This is the question of when a reaction 
can be added to or removed from a mechanism without changing the computed perform- 
ance significantly. An attempt to answer this question has been made by comparing the 
reactions in a given system with each other. To do this we use the following ratio for 
each reaction: 



( 6 ) 
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Here (x^ ^ is the largest Xtj . value for all the reactions in the mechanism. 

V H ’J'max H ’ 3 

The smaller this ratio is for a given reaction, the less important it should be in deter- 
mining performance. There ought to then be a limiting value of this ratio that has the fol- 
lowing property: all reactions having Ry j values below this limiting value will change 
kinetic specific impulse by no more than ±0. 05 (lb)(sec)/lb (4. 9 (N)(sec)/kg) when added to 
or removed from a reaction mechanism. The three sets of Xy j values given in 
table IV were used to compute R w a values for the three assumed conditions, and the 
results are shown in table X. Additional performance calculations were then made for 
these three cases and a few others not shown. Various reactions in the mechanism were * 

systematically eliminated, starting with the least important ones, in an attempt to get an 
approximate value of the limiting Rtt • ratio. These computations gave numbers be- * 

tween 0. 015 and 0. 025 as the limiting R R . values. These few calculations cannot be 
used to establish a good quantitative criterion for identifying unimportant reactions. The 
results do give, however, an approximate guide to reaction importance. Any reaction 
with an R H • value greater than about 0. 02 should probably be included in a mechanism. 
However, reactions with Ry j less than about 0. 02 may probably be eliminated without 
causing more than a ±0. 5 (lb)(sec)/lb (4. 9 (N) (sec) /kg) change in the computed kinetic 
specific impulse. It was found previously (table V) that the reaction F + Hg ~ HF + H 
could be dropped from consideration at low O/F but not at higher O/F values. This re- 
sult is consistent with the observed change of Rtt * with O/F shown in table X for this 

J 

reaction. We may also determine whether Ry j values are consistent with the obser- 
vation that the addition of reaction (XIV) (F + OH s HF + O) has no significant effect on 
computed performance. The information for cases 1 and 3 in tables VHI and IX gives the 
following Ry j values for reaction (XIV) and corresponding impulse changes AI yac 
when this reaction is added to the original mechanism. For case 1, Ry ^ = 0. 002 and 
AI = 0. 23 (lb)(sec)/lb(2. 3 (N)(sec)/kg), and, for case 3, R H • = 0. 14 and AI = 0. 51 
(lb)(sec)/lb(5. 0 (N)(sec)/kg) . Thus, case 1 gives alow AI vac consistent with the Ry j 
value much lower than 0. 02; case 3, however, gives a AI vac smaller than would be ex- 
pected for an Rtt • value much larger than 0. 02. It appears that the stated rule for re- 
action importance will indicate the inclusion of a reaction whose effect on performance is 
very close to our arbitrary limit for significant effect on performance. These computa- 
tions have shown that reaction (XIV) might have to be considered under some assumed 
conditions, but not those used in the present work. 

Comparison of experimenta l and theoretical results. - The experimental measure- 
ments reported in reference 6 were compared with the present theoretical nonequilibrium 
calculations after correction by the method described in the Procedure section. Two of 
the 11 test results given were not used, because they were stated to have error in the 
thrust measurement. The nine data points used are the same ones used for a similar 
comparison in reference 7 and shown in figure 123 of that work. The typical correction 
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TABLE X. - REACTION RATING PARAMETER RATIO 8, VALUES FOR METHANE 
FLUORINE - OXYGEN NOZZLE RECOMBINATION REACTIONS 


[Chamber pressure, 100 psia (6. 895X10 5 N/m 3 ) . ] 


Reaction 

number 

Reaction 


3 

Oxidant-fuel mixture ratio, O/F 




6 

6 

Fluorine content in oxidant, wt. 

% 





82.6 

73. 5 

82.6 



Reaction 

Reaction 

Rank 

Reaction 

Reaction 

Rank 

Reaction 

Reaction 

Rank 



rating 

rating 


rating 

rating 


rating 

rating 




parameter, 

parameter 


parameter, 

parameter 


parameter, 

parameter 




X H,j 

ratio, 


X H, j 

ratio, 


x h, j 

ratio, 





r hj 



R H,j 



P H,j 


I 

H+F+MitHF+M 

-7 627 

0. 0766 

2 

-124 600 

1.000 

1 

-232 600 

1.000 

1 

n 

H+H+Mrfl 2 +M 

-99 610 

1. 000 

1 

-20 590 

. 165 

4 

-24 850 

. 107 

2 

m 

F+H 0 itHF+H 

-731 

. 00734 

3 

-24 600 

. 197 

3 

-18 350 

.0789 

3 

IV 

CH+CH^C 2 H 2 

-2. lXlO" 3 

2. 1X10" 8 

5 



— 



— 

V 

C+H+JVfcCH+M 

-1.36 

1.37x10 0 

4 



-- 



— 

VI 

cch-o+m^co 2 +m 



— 

-699 

.0056 

10 

-74 

.00032 

8 

vn 

COfOH^C0 2 +H 



-- 

-5 000 

.0401 

7 

-353 

.0015 

7 

vm 

H+OH+M^HgOfM 



— 

-56 640 

.455 

2 

-3 946 

.0170 

4 

IX 

o+o+m^o 2 +m 



-- 

-2 951 

.0237 

8 

-31 

.00013 

9 

X 

OfH+M^OH+M 

j 


-- 

-9 129 

.0733 

5 

-1 004 

.0043 

5 

XI 

H 2 +OH^H 2 CM-H 



-- 

+6 675 

. 0536 

6 

+511 

.0022 

6 

XII 

h+o 2 ^oh+o 



-- 

-2 414 

. 0194 

9 

-24 

. 00010 

10 

xm 

oh 2 ^oh+h 

1 



-- 

+ 194 

.0016 

11 

-18 

. 00008 

11 
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from the original experimental I vac to the kinetically limited value is about 40(lb)(sec)/lb 
(392 (N) (sec) /kg), in the present work. Of this total, the correction for combustion effi- 
ciency is approximately 35 (lb)(sec)/lb(lb(343 (N)(sec)/kg) and the remainder is the cor- 
rection for heat transfer, friction, and nozzle divergence losses. The experimentally de- 
rived kinetic performance data are plotted as the circled points in figures 5(b) and 6(b). 

In figure 5(b) the computed curve using the reference rate constants gives the best agree- 
ment with the data points . However, in figure 6(b) a quantitative calculation had to be 
done to determine whether the data points are fitted better by the reference curve or the 
curve using the reduced rate constant for the H + F reaction. To measure the deviation 
from each of these two curves, the percentage standard deviation of the data points was 
calculated using the formula 



where I is the corrected experimental value and n is the number of data points. The 
value of a for the curve using reference rates is 1. 7 percent, while a is 1.4 percent 
for the curve with the reduced value of kj. This means that either set of rate constants 
gives a theoretical curve that matches the corrected experimental results to about the 
same accuracy. Therefore, using the reference rate constants presented herein has 
given a theoretical performance curve that agrees with a particular set of corrected ex- 
perimental data to a percentage standard deviation of about 1-^ percent. The closeness of 
the agreement may be fortuitous, in view of experimental uncertainties as well as uncer- 
tainties in the method used to correct the experimental data. This method is simpler and 
less exact than the one used in reference 7; that work compares theoretical and cor- 
rected results which are both different from the results presented herein. 


Butene-l-Fluorme - Oxygen System 

The mechanisms and reference reaction rate constants given previously have been 
used to calculate the kinetic performance of the butene -1 -fluorine - oxygen propellant 
system. The results are presented for a chamber pressure of 100 psia (6. 895x10^ 
N/m^) and an O/F range of 2. 5 to 5. For this combination, oxidant mixtures ranging 
from 54. 3 to 82. 6 weight percent fluorine were used. An oxidant having 70. 4 weight 
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400 


k 



2.4 2.8 3.2 3.6 4.0 4.4 4.8 5.2 

Oxidant-fuel ratio, O/F 


(a) Fluorine content in oxidant, 54.3 weight percent, stoichio- 
metric O/F, 4.99. 



(b) Fluorine content in oxidant, 70. 4 weight percent, stoichio- (c) Fluorine content in oxidant, 82. 6 weight percent; stoichio- 
metric O/F, 5.76. metric O/F, 6.56. 

Figure 7. - Butene- 1-fluorine - oxygen performance. Chamber pressure, 100 psia (6.895X10 5 N/m); area ratio, 40. 


percent fluorine corresponds to the chemical equation 

C 4 Hg + 20 2 + 4F 2 - 4CO + 8HF (B) 

which should give maximum performance. Figure 7 shows complete theoretical results 
(equilibrium, frozen, and kinetic) for three different oxidizers. In figure 8 the equili- 
brium and kinetic data are cross plotted to show the effect of oxidizer composition on 
performance. These curves show that the maximum kinetic performance of 385 pound- 
seconds per pound is obtained at an O/F value of 3 . 6 for an oxidizer containing about 
70 weight percent fluorine, in agreement with reaction (B) . A few experimental data 
points are reported in reference 6. These were corrected as described previously and 
are plotted in figure 7(b). Three of the points are in very good agreement with the com- 
puted curve, and the standard deviation of all the data points is 0. 8, or about 1 percent. 

A comparison of the performance of the methane fluorine - oxygen and butene- 1- 
fluorine - oxygen systems is shown in figure 9. Kinetic and equilibrium curves are 
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Figure 8. - Effect of oxidant composition on performance of butene-1- 
fluorine - oxygen. Chamber pressure, 100 psia (6.895xl0 5 N/m 2 ); 
final area ratio, 40. 
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Figure 9. - Comparison of performance of methane-fluorine - oxygen and butene- 
1-fluorine - oxygen propellant systems. Chamber pressure, 100 psia 
(6.895X10 5 N/m2) ; area ratio, 40. 



plotted using the oxidizer that gives best performance for each system. The maximum 
kinetic specific impulse for butene-1 is about 385 (lb)(sec)/lb (3776 (N)(sec)/kg), where- 
as, for methane, it is 398 (lb)(sec)/lb (3903 (N) (sec) /kg). This comparison on the basis 
of specific impulse shows that methane is definitely superior to butene- 1 when combined 
with the fluorine -oxygen oxidizer. 


CONCLUDING REMARKS 

Two general conclusions can be drawn from the results of the present work. First 
of all, three, or at most four, chemical reactions are really rate controlling in the over- 
all recombination process of the multicomponent gas mixtures studied. Although a dozen 
or more individual reactions may actually occur, only the rate constants of a few three - 
body recombination processes have to be known to better than one order of magnitude un- 
certainty for computing kinetic rocket performance within an uncertainty of ±0. 5 pound- 
seconds per pound (4. 9 (N)(sec)/kg). Many other reactions must be present in the mech- 
anism, but their rate constants can be changed in most cases by an order of magnitude or 
more without changing the computed specific impulse more than this uncertainty. The 
role of these reactions, which are the bimolecular processes, is primarily to maintain 
concentrations of the free radical and atomic species that combine in the rate- controlling 
reactions. 

In the hydrocarbon-fluorine - oxygen systems studied herein, the rate controlling 
reactions are H + H + Mj Hg +M, H+F+Mi HF + M, and, under a few conditions, 

H + OH + M s HgO + M. The bimolecular reaction occurring to the greatest extent is 
F + Hg s HF + H. However, performance is not significantly changed when its rate con- 
stant is changed from 1/10 to 10 times the reference value. 

The other significant idea that emerges from this work is verification of the ability 
of a one-dimensional idealized flow model to calculate performance in good agreement 
with experiment when a selected set of reaction rate constants is used. This agreement 
was obtained in reference 2 for the hydrogen-fluorine system. In the present work, com- 
puted kinetic performance matches the corrected experimental results with a standard 
deviation of about l^- percent for the methane-fluorine - oxygen and 1 percent for the 
butene- 1 -fluorine - oxygen system. One should not conclude, however, that the one- 
dimensional theory will necessarily describe other systems equally well, because uncer- 
tainties still exist about the exact manner of correcting the experimental measurements. 


SUMMARY OF RESULTS 


Nonequilibrium rocket performance calculations have been performed for the 
methane -fluorine - oxygen and butene- 1 -fluorine - oxygen systems for chamber pres- 
sures from 100 to 300 psia (6. 895X10 5 to 2. 069X10 6 N/m 2 ) and oxidant-fuel weight ratios 
from 2. 5 to 6. A truncated perfect nozzle profile with an area ratio of 40 was used. 
Computed equilibrium and nonequilibrium specific impulse are reported and have been 
compared with experimentally derived kinetic performance data. Two different rating 
parameters have been used to quantitatively rank all chemical reactions in order of their 
importance in controlling performance. The predictions of these parameters were 
checked by individual rate -variation calculations. The results of this investigation are 
as follows: 

1. Theoretically computed performance using a one -dimensional theoretical model 
and a selected set of reaction rate constants agrees with "corrected” experimental data 
with a standard deviation of about percent for methane and 1 percent for butene- 1 
combined with fluorine -oxygen mixtures. 

2. Two rating parameters (X. proportional to net reaction rate and X H • propor- 

J J 

tional to net heat exchange rate), predict that the same four reactions are most signifi- 
cant for the conditions used. These are H+H + M; Hg +M, H + F+ Ms HF + M, 

H + OH + M 2 H 2 O + M, and the bimolecular reaction F + Hg £ HF + H. The parameter 

Xtt • is somewhat better than X. for quantitatively ranking the important reactions on 

n,] J 

the basis of performance change caused by uncertainty in each reaction rate constant. 

3. The maximum theoretical nonequilibrium performance for methane -fluorine - 
oxygen is 398 pound-seconds per pound (3902 (N)(sec)/kg) at 100 psia (6. 895x10^ N/m 2 ) 
chamber pressure and an oxidant-fuel ratio O/F of 5. 1 with the 82. 6 -weight-percent 
fluorine oxidizer. For butene -1 -fluorine - oxygen, the maximum performance at 100 psia 
(6. 895X10 5 N/m 2 ) and O/F = 3.6 is 385 pound-seconds per pound (3776 (N)(sec)/kg) 
using the 70. 4 -weight-percent oxidizer. 

4. Calculated performance is most sensitive to the H + H + M s Hg + M reaction rate 
at low O/F values for methane-fluorine - oxygen using both the 73. 5 and 82. 6 weight 
percent fluorine oxidizers. At high O/F values, the reactions H+F+M; HF + M and 
H + OH + M s H 2 0 + M are both rate determining for the 73. 5 weight percent fluorine 
oxidizer. However, increasing the fluorine content to 82. 6 weight percent at high O/F 
eliminates the water recombination as an important reaction in the mechanism. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, September 19, 1968, 

129-01-02-01-22. 
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APPENDIX A 


SYMBOLS 



eq 

X fr 


vac 

k-i • 

b,J 


k f,i 


M 

O/F 

(°/f) st 

P c 

R 




pre -exponential factor in rate constant expression 

3 3 

concentration of species i, (lb)(mole)/ft ; (g)(mole)/cm 

activation energy, cal/mole 

impulse recovery factor defined by eq. (3) 

heat of reaction at 298 K, kcal/mole; kj/mole 

experimentally derived, kinetically limited vacuum specific impulse de- 
fined by eq. (2), (lb)(sec)/lb; ((N)(sec)Ag) 

theoretical equilibrium vacuum specific impulse, (lb) (sec) /lb; (N)(sec) /kg 
theoretical frozen vacuum specific impulse, (lb)(sec)/lb; (N)(sec)/kg 

theoretical kinetically limited vacuum specific impulse, (lb) (sec) /lb; 

(N) (sec) /kg 

theoretical vacuum specific impulse, (lb)(sec)/lb; (N)(sec)/kg 

O 

reverse reaction rate constant for j reaction, cm°/(mole)(sec) or 
cm®/ (mole) ^(sec) 

th 3 

forward reaction rate constant for the j reaction, cm /(mole) (sec) or 
cm®/ (mole) ^(sec) 

third-body molecule in recombination reaction 
assigned oxidant to fuel weight ratio 
stoichiometric O/F value 

2 

chamber pressure, psia; N/m 

assigned equivalence ratio = (O/F) ./(O/F) 


reaction rating parameter ratio = 


X, 




( X H.i)max 


reaction rating parameter based on net enthalpy exchange rate for j 
action defined by eq. (5), (kcal) (ft®) / (lb) (sec) ; (kg) (m®) /(kg) (sec) 

largest X„ . value for all reactions in a given mechanism 

ii, ] 


.th 


re- 
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1 


X. 

J 


1?C* 


"ij 

1i 


reaction rating parameter based on net reaction rate for j 
eq. (4), (mole) (ft 3 ) /(lb) (sec); (mole) (m 3 ) /(kg) (sec) 

volumetric conversion rate for reaction, mole /(ft 3 ) (sec); 

characteristic velocity efficiency 

stoichiometric reaction coefficient for j chemical reaction 

fVi 

stoichiometric product coefficient for j chemical reaction 

3 3 

gas density, lb/ft ; g/cm 


Subscripts: 
corr corrected 
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reaction defined by 

O 

mole/(m )(sec) 


J 



APPENDIX B 


CALCULATION OF CONDITIONS FOR MAXIMUM THEORETICAL PERFORMANCE 
OF HYDROCARBON-FLUORINE - OXYGEN PROPELLANT 

The propellant formula is written using the notation of reference 9 

CH y + ^(OF^) 

Then the oxidant-fuel ratio of this mixture is given by 

O _ <*(16 +19/3) ( B1 

F 12.011 + 1.008 y 

Now, when it is assumed that the propellant composition is assigned, /3 and y are 
known. We consider two cases for finding a. 

(1) Fluorine -oxygen atom ratio less than or equal to hydrogen- carbon atom ratio 
(/3 < y): The only products are assumed to be HF, CO, and HgO. The hydrogen atom 
balance then requires that 


a = 21±1 (B2) 

P+2 

Thus, the complete propellant composition is known, and the oxidant-fuel ratio for maxi- 
mum performance using a particular oxidant is 

ZO\ = (y + 2) (16 + 19 /3) (b3) 

\F/ max (J8+ 2)(12. 011+ 1.008 y) 

Notice that <*, the oxygen-carbon atom ratio, is greater than 1 and that this case 
corresponds to reaction mechanism B of this report. 

(2) Fluorine -oxygen atom ratio greater than hydrogen- carbon atom ratio (P > y): The 
final products are assumed to be carbon (C), HF, and CO. For this reaction the hydro- 
gen and fluorine atom balances require that 


a = ?- 




(B4) 
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The complete propellant composition is now known and the oxidant-fuel ratio for maximum 
performance is given by 


/ox = r(16 + 19 /3) 

VF/ max /3(12.011 + 1.008 y) 

For this case the oxygen-carbon atom ratio is less than 1 and this corresponds to reaction 
mechanism A. 

It should be mentioned that equations (B3) and (B5) give (O/F) values slightly 

max 

different from those actually obtained from equilibrium computations. This is because 
small amounts of products other than those assumed are present in the gas mixture. 
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